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Abstract 

Using a displacement operator formalism, I analyse the depopulation of the 
vibrational ground state of trapped ions. Two heating times, one characterizing 
short time behaviour, the other long time behaviour are found. The short time 
behaviour is analyzed both for single and multiple ions, and a formula for the 
relative heating rates of different modes is derived. The possibility of correction 
of heating via the quantum Zeno effect, and the exploitation of the suppression of 
heating of higher modes to reduce errors in quantum computation is considered. 
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Individual or multiple ions can be confined in a radio-frequency Paul trap and using 
sophisticated laser techniques cooled to the quantum mechanical ground state [jl], |2|. 
Such systems allow experimental preparation and measurement of non-classical motion 
states of the ions, and are therefore of great current interest in physics [[|. Further- 
more important technological applications for such systems, such as the practical im- 
plementation of quantum computation, have recently attracted considerable attention 
[!> H) H) 01 • Quantum information can be stored in the internal quantum states of the 
ions (which constitute the quantum bits, or "qubits" of the computer), and, using ultra 
narrow bandwidth lasers, quantum gate operations can be realized between pairs of 
qubits using quantum states of the collective motion of the ions in the harmonic con- 
fining potential as a quantum information bus. If this bus were to become degraded by 
heating information would be lost, and so it is of great importance to maintain the ions 
in their motional ground state as long as possible. Of the many practical roadblocks 
standing in the way of success in realization of an ion trap quantum computer, one 
of the most important is the very fragile nature of this motional ground state, due to 
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interactions with various ambient electromagnetic fields. In this letter I present a the- 
oretical analysis of the depopulation of the motional ground state of ions due to such 
fields, which can be loosely characterized as heating (although it should be stressed 
that in this analysis relaxation to a thermal distribution is not considered). 

Various analyses of decoherence mechanisms in ion trap quantum computers have 
been carried out || 0, [TO, O, O. Such effects as spontaneous emission from the 



internal states || [10], pjl , or dephasing of the internal states due to the ions' zero-point 
motion ||, both of which effects can degrade or destroy quantum information, have 
been considered. Also the relaxation of an ion to a thermal state has been considered 



using a perturbation method, valid for describing long time behavior other work 
has been done on the effects of laser amplitude and phase stability |L4|, [15] . 

In this letter I adopt a rather different approach than that normally used in analysis 
of decoherence problems. When the electromagnetic fields causing the heating can be 
treated classically, the equations of motion reduce to a solvable system, namely a 
simple harmonic oscillator driven by a classical field. One can thus dispense with 
the density matrix/master equation formalism, using instead the interaction picture 
to solve exactly for the wavefunction describing the state of the ion. One can thus 
obtain an expression for the population remaining in the ground state as a function of 
time. An average of this population over an ensemble of realizations of the classical 
perturbing field (assumed to have Gaussian statistics) yields simple expressions for the 
time dependent behavior of the depopulation both for short and for long time limits. 

In the interaction picture the Hamiltonian describing the interaction of a single 
ion of mass M with a classical electric field E(t) is given by the following formula 

H = -eE(t)x 



ih 



u(t)tf -u*(t)a] (1) 



where x is the operator for the position of the ion and u(t) = ieE(t) exp(ioj t) / \ / 2MTiuj , 
e being the ion charge and % Planck's constant divided by 27r. I am considering only 
the motion of the ion along the weak confinement axis of an anisotropic trap, of the sort 
suitable for quantum computation; thus E(t) represents the component of the electric 
field along that axis. The harmonic binding potential is characterized by the angular 
frequency u . The operators a (a*) are the zero time annihilation (creation) operators 
for the harmonic motion of the ion in the harmonic well. 

The dynamics of such a driven quantum harmonic oscillator can be solved exactly 



HI, O] . The wave function at some instant t is related to the initial wave function by 



the following expression 

\il>(t))=exp[i<l>(t)]D [v(t)\ h/>(0)>, (2) 

where D[v] = exp va) — v*a is the displacement operator for coherent states of a 
harmonic oscillator, (f)(t) is a phase factor (which turns out to be unimportant for the 
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current problem) and the amplitude v(t) is given by the formula 

v(t) = f u(t')dt' 
Jo 



y/2Mhu Jo 



E(t')exp(iu t')dt'. 



(3) 



As a figure of merit, let us introduce an average fidelity of the ground state, defined 
by the following formula: 

F(t) = (\(m\mu 2 )f, (4) 

where (. . .) q denotes a quantum mechanical average, and (...)/ an average over an 
ensemble of realizations of the classical random field E(t). If one assumes that the ion 
is initially in the ground state, then its state evolves into a coherent state of amplitude 
v(t). Thus the probability of remaining in the ground state can be found in closed 
form, and the fidelity of the ground state is given by the formula 



F(f) = (exp -Kt)| 2 }/. 



(5) 



If one assumes Gaussian statistics for the classical random field E(t) |l8j (which can be 
justified by assuming that the field is due to many uncorrelated random sources, and 
then invoking the central limit theorem), the average over the field ensemble can be 
determined by performing an integration using the appropriate probability distribution. 
The result is as follows: 



F(t) 



l + 2{\v{t)\% 



-1/2 



(6) 



+ {(Hm f ) -\(v(tHt)} f 

An alternative measure of the heating is the mean excitation number, defined by 

n(t) = ([m)\^^\m),])f- (7) 

Using eq.fl^) and eq.(11.3-13) of ref. HT7| , this can be rewritten as 

n = ([(0\n + v*(t)a + tfv(t) + \v(t)\ 2 \0} q }} f 

= (Ht)\ 2 )f. (8) 

The correlation functions appearing in eqs.(^) and (|8[) can be evaluated using 
eq.(|3|). Using the symmetry property of the autocorrelation function, one obtains: 



<Kor>/ 



2.2 



x)jE(xt) cos (xu t) dx 



X 



Q 2 t 
i 

7£i(xt) sin [(1 — x)u) t] dx 



(9) 



(10) 



D □ 



Figure 1: The fidelity of the ground state as a function of time, illustrating the results 
given in eqs.(^), ( PD and ([12]). The parameters used were as follows: curve (a) u T=l, 
u)qTi= 1; curve (b) u T=l, u)qTi= 8.5; curve (c) u Q T=l, u Ti= 41; curve (d) u T=l, 
Uqt^ 128.5. 
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Two limiting cases are of interest, namely the behavior for short times and for 
very long times. For short times (i.e. i<T, l/^o) the mean excitation number and 
fidelity are given by 



n(t) pa 1 



F(t) 



( 




t 2 + O (t 3 ) 



(14) 



This result, i.e. that for short times the "decay" of the ground state population is non- 
exponential, allows the possibility of maintaining the ion in its ground state via the 
quantum Zeno effect |L9|, |20| . If repeated measurements of the ground state population 
can be made on time scales much shorter than 1/uq, then the ion will in principle 
remain in the ground state for a much longer time. However, as the quadratic behaviour 
persists for a short time only (see Fig.l), such a technique would be very difficult to 
implement. Furthermore, by invoking the time-energy uncertainty principle, one can 
see that measurements carried out on time scales much less than 1/uq will result in 
excitation of higher number states. Hence it is unlikely that the Zeno effect will be a 
practical method of error correction in these circumstances. 

One can find asymptotic expressions for F(t) and n(t)m the long time limit t ^> 



where t = T(l — UqT 2 )/(1 + uj^T 2 ). These results are in qualitative agreement with 
that obtained earlier by Lamoreaux |H| using a perturbative density matrix approach. 

In order to gain some insight into the order of magnitude of the heating effect, 
one can calculate a formula for the value of T\ when the ambient electric field E(t) is 
thermal, assuming a correlation time T pa h/KsQ <C l/^0) (where Kb is Boltzmann's 
constant and 9 is the temperature): 



where a is the Stefan-Boltzmann constant and eo is the permittivity of free space. 
Using data from the mercury ion experiment described in ref.fllj (i.e. M = 3.29 x 1CT 25 
kg, uq = (27r)4.66MHz), a heating time Ti=135 ms implies an effective temperature of 
T— 4.6 K. This result should be treated with caution: the assumption that the random 
ambient field in the vicinity of the ion is thermal is not justified. 

Let us now consider the case of multiple ions confined in a linear configuration. 
Because the ions are interacting via the Coulomb force, their motion will be strongly 
coupled. The small amplitude fluctuations are best described in terms of normal modes, 
each of which can be treated as an independent harmonic oscillator [0. If there are N 
ions in the trap, there will be a total of N such modes. I will number these modes 



T, l/u : 




(15) 



(16) 



?>C6qUJqMKb 



(17) 
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in order of increasing resonance frequency, the lowest (p = 1) mode being the center 
of mass mode, in which the ions oscillate as if rigidly clamped together. In the quan- 
tum mechanical description, each mode is characterized by creation and annihilation 
operators al and a p (where p — 1, . . . N). The Hamiltonian in this case is given by the 



expression 



N 



H 



n=l 
N 



ih ^2 [u P (t)al - u*(t)Of 



18) 



P =i 



where E n (t) is the electromagnetic field at the n-th ion of the string of ions, x n (t) is 
the operator for the displacement of the n — th ion from its equilibrium position and 



u p {t) 



is 



N 

Entyb^ exp (iy/wuot) . 



(19) 



In eq. ( |19D , is the n-th element of the p-th normalized eigenvector of the ion coupling 
matrix [0, u p being its eigenvalue. Again the evolution of the state of the ions can be 
solved exactly: 



N 



= exp [z$(*)] n D p [v p (t)} |^(0)), 
P =i 



(20) 



where D p [v p } = exp 



VpClp Vp(i p 



is the displacement operator for the p-th mode, and 



v p (t) = / u p (t')dt'. 



(21) 



As before, one can find a closed form expression for the fidelity of the ground state of 
the string of ions. For multiple ions, the mean excitation number is given by a formula 
analogous to eq.(16), where the characteristic decay time is given by: 



TN = Ti 



N N Up)Up) 

p=l m,n=l V 



-1 



(22) 



Imn being the degree of coherence of the field at the positions of ions n and m evaluated 
for zero time delay |17j (the field having been assumed to be cross-spectrally pure 



21], P2[), and the coherence time T has been assumed to be much less than 1/uq- 
In the coherent limit (7 mn = 1), the formula for reduces to the following simple 
expression 

tn = p (23) 
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Figure 2: The heating time of the ground state for different numbers of ions, in the 
spatially coherent (plain curve) and spatially incoherent (dashed curve) limits for the 
ambient electric field. 



One can also consider the mean excitation numbers of the different modes. These 
quantities are given by expressions analogous to eq.(|8|), with the heating time of the 
p-th mode, when there are N ions in the string, given by 



r N ,p = N — n. (25) 

/ , u m u n Imn 
m,n=l 

In the coherent and incoherent limits, the heating times of the different modes are 
given 

tnv = I ^ \ coherent (26) 

' p [ oo p > 1. v ' 

t n, p = y/J^pTi incoherent. (27) 
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The heating times for the coherent case is potentially a very important result. Only 
the lowest (p — 1) centre of mass mode will be heated up by spatially coherent fields. 
Thus the state of the ion oscillations is given by: 

= exp#(t) <g> |0) a ® |0) 3 ® . . . ® |0>jv} , (28) 

the modes other than the center of mass mode remaining in their ground states. It 
should therefore be possible to perform quantum logic operations without degradation 
due to heating, simply by utilizing these higher modes as the quantum information bus 



rather than the center of mass mode [23 
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